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Central  Limit  Theorem  for 


Parking  Models  on  the  Line  and  Plane 
by  Howard  J.  Weiner 

University  of  California  at  Davis  and  Stanford  University* 

1.  Intr oduc t i on .  Three  models  for  random  parking  (=  packing)  of  line 

segments  on  a  longer  line  and  correspondingly,  squares  on  a  larger 
rectangle  are  given,  where  for  each  model,  either  the  Renyi  or  the 
Solomon  parking  method  is  used. 

The  one -dimensional  first  moment  results  are  given  in  Solomon  (1966), 
(1970),  Weiner  (1978). 

In  the  alternating  car  size  model,  a  car  of  length  a  is  placed 
uniformly  at  random  along  a  curb  (segment)  of  length  x.  It  is  considered 
praked  there.  A  second  car,  of  length  b  is  placed  uniformly  at  random 
along  Che  curb  and  parked  if  and  only  if  it  does  not  overlap  the  a-car. 

This  is  the  Renyi  parking  mechanism.  In  the  Solomon  parking  mechanism, 
a  newly-placed  car  which  overlaps  an  already-parked  car  is  moved  the 
shorter  distance  so  as  to  be  immediately  adjacent  to  the  overlapped  car 
on  the  curb,  and  is  parked  there  if  and  only  if  it  does  not  now  overlap 
yet  another  parked  car,  or  the  boundary  of  the  segment.  If  a  b-car  is  not 
parked,  it  Is  discarded,  and  another  b-car  is  placed  I.I.D.  at  random  as 
the  other  b-car.  If  it  is  parked,  then  an  a-car  is  placed,  and  the  process 
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continues  until  no  further  cars  can  be  parked.  If  the  b-car  is  not 
parked,  again  it  is  discarded  and  another  b-car  is  placed,  and  again 
the  process  continues  until  no  more  cars  may  be  parked.  This  process 
may  be  similarly  defined  for  the  plane  for  aXa  and  bxb  squares, 
respectively,  with  either  of  the  Renyi  and  Solomon  mechanisms. 

The  abacus  model  consists  of  a  rectangle,  (k+l)xx,  of  k  horizontal 
lines,  each  one  unit  apart.  The  width  of  the  rectangle  (the  length  of 
each  line)  is  x.  Assume  a,  b  are  integers.  A  line  is  chosen  uniformly 
at  random  from  the  integers  {l,2,...»k],  and  centered  on  this  line,  an 
axa  (for  definiteness)  square  Is  placed.  Next,  a  bxb  square  is  selected 
and  placed  uniformly  on  a  line,  and  parked  if  there  is  no  overlap.  If 
there  is  overlap  with  the  axa  parked  car  on  a  given  line  (Renyi  model) 
or  if  a  car  to  be  placed  is  vertically  adjacent  to  some  portion  of  a  car, 
without  a  line  of  space  between  cars,  the  bxb  car  is  placed  I.I.D.  as 
before  until  It  is  parked.  Then  another  axa  car  is  placed  I.I.D.  and  the 
process  continues  until  no  further  cars  may  be  parked.  The  Solomon  model 
allows  that  if  two  cars  overlap  on  the  same  line,  the  one  to  be  placed  may 
be  moved  adjacent  horizontally  (minimal  motion)  to  the  already-parked  car 
and  parked  if  there  is  now  no  overlap  with  another  parked  car,  if  the 
boundary  is  not  overlapped,  and  if  at  least  one  line  of  vertical  space 
exists  between  the  car  and  already-parked  cars  on  other  lines. 

The  random  car  size  model  to  be  considered  has  DXD  square  cars  to 
be  parked,  either  on  a  rectangle  or  abacus  with  the  Renyi  mechanism,  where 
D  is  chosen  from  a  distribution  F,  with  density  f  (x)  >  0,  all  x  >  0. 

If  a  DXD  car  fits  the  boundaries  of  the  given  rectangle  and  does  not  overlap 
another  already-parked  car,  (and  leaves  enough  vertical  space  in  the  abacus 
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grid  case),  it  is  parked,  otherwise  it  is  discarded  and  the  process  stops. 
If  the  car  is  parked,  another  D  is  chosen  I.I.D.  from  F  and  a  DXD  car 
is  placed  I.I.D.  as  the  other  cars,  and  if  It  is  not  parked  as  described 
above,  the  process  stops. 


2.  One-Dimensional  Models.  Define,  for  a,  b  «  x, 

(2.1)  X^(x)  (X^(x))  =  the  random  variable  of  total  number  of  cars 

which  may  be  parked  on  a  curb  of  length  x  in  the  alternating 
a  and  b  length  car  Solomon  model,  starting  with  an  a-car  . 
(b-car) 

K  (x)  £  E(X  (x)) 

a  a 

K^x)  h  ECX^x)). 

By  considerations  as  in  (Solomon  (1966)),  it  may  be  shown  that 


(2.2) 


-lim_X'1(^(x))  ~  C 


Pl-\ 

P2^ 


P- 


Define 


(2.3) 


VO 

L  (s)  «=.  f  e  SXK  (x)  dx 


a  '  J  0  a 


^(b)  »  J^e'SXKb(x)  dx 


that  for  some  constants  0^,  02,  6^  6^  which  can  be  determined,  for  s  i  0 
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(2,4) 


Lfl(s)  ~  Pj^s  2  +  Sj^b"1 

-2  -1 

LjjC8)  ~  P28  +  ^28 

Let 

(2.5)  Xa(x)  s  Pj*  +  6i 

^b(x)  s  P2x  +  f>2  * 

It  may  be  shown  that  £(x)  satisfies  the  same  integral  equation  as 
_  f  Ka^X^  '\ 

£(x)  =  V  K^x)  )'  HenCe 

(2.6)  Xa(x+a)  -  Xa(x+a)  -  (^(*) ”\ (x)  )JC[  (_a> 0)|j(x,x+a) ] 


+  [(X^cty-J^ct))  +  (Xj;(x-t)-£b(x-t))] 

and  a  similar  equation  for-X^x+b)  -  X^(x+b),  where  the  a,  b  are  inter¬ 
changed,  where  t  is  chosen  uniformly  on  [0,x],  X'(t),  X"(t)  are  I.I.D. 
as  X  (t),  conditional  on  t.  Squaring  both  sides  of  (2.6)  and  taking 
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(2)  2 

expectations,  with  IT  (x)  E  E(X  (x)-jt  (x))  , 

O  A  A 


(2.7) 


£**><«>  - 


+  sk  dt* 


and 
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M^2)(x+b)  =  -~r  M(2)(x)  +  — ~  f  M(2)(t)  dt  +  — | ~  f  m  (t)m  (x- 
b  v  '  x+2b  a  v  x+2b  j  _  a  x+2b  J  _  a  a 


where 

(2.8) 


m  (t)  s  E(X  (t)-£  (t)), 

a  a  d 


Define 


“b 


(t)  =  e (x,  (t)  -  je  (t)> 


(2.9) 


La(S)  h  j  e~SXM^2)(x)  dx 


dx 


r  (s)  s  f  e_SXm  (t)  dx 
a  a 

a 


CD 

rb(s)  s  J  e”sbmb(x)  dx 


L  = 


La(s) 

Lb(s) 


Taking  Laplace  transforms  in  (2.7),  one  obtains 


(2.10) 


L'  +  A  L  =  - B 


where 


t)  dt. 
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B  = 


Define 


(2.11)  K^a(x) (K  (x))  =  mean  total  number  of  a-length  (b-length) 
cars  which  can  be  parked  on  [0,x) . 

Assume  a  <  b. 

Clearly 

(2.12)  Klb(x)  <  Ka(x)>Kb(x)  <  Kla(x), 

and  a  similar  inequality  holds  for  the  respective  non-normalized  second 
moments . 

Hence  we  seek  a  solution,  either  by  power  series,  or  Volterra-like 
integral  equation  representation  to  (2.10)  such  that 


(2.13) 


lim  s  (L.(s)  =  a.  >  0,  i=l,2. 

sjO  1  1 


(2)  (2) 

Hence  Ma  (x) ,M^  (x)  are  respectively  expressible  as  a  sum  of  three 

terms,  each  of  which  is  increasing  in  absolute  value.  By  Abelian,  Tauberian 
theorems  (Widder  (1946),  p.  182  and  p.  192), 


(2.14) 


-ir“a2>(x)V  ..  -lAVar  Xa(x>^ 

lim  x  .  .  J  =  lim  x  (  !  = 

x V  Mb  (x)  x ->  co  Var  X^x)-^  O' 


O) 
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For  the  abacus,  define 


(2.15)  Y  (k)  =  total  number  of  alternately-placed  a  and  b  segments 
a 

on  a  k-line  one-dimensional  abacus  starting  with  an  a 
segment. 

Y^(k)  is  the  similar  quantity  starting  with  a  b  segment. 


One  may  conclude  by  similar  considerations  as  in  (Solomon  (1966)) 


that 


(2.16) 


Yb(k)y 


X(k) 


for  large  k,  where  X  (k) ,  JL  (k)  are  linear  functions  of  k. 

a  d 

Denoting 


(2.17)  Mf2)<k)  =  E(Y  (k)-X  (k))2 

a  d  d 

M^2)(k)  =  E(Yb(k)-Xb(k))2 
in  (k)  =  E  (Y  (k)  -X  (k)  ) 

Pi  a  a 

^(k)  =  E(Yfa  (k)-Xb(k)  ) 


it  follows 

(2.18) 


that  with 

K(2)(k+a) 

a 

H^2)(k+b) 


k  an  integer,  k  »  a  or  b. 

2  k  2  k  (2) 

-  i  -  k  \o 

2  k  2  k  (2) 

-r  S  m  (X)m  (k-X)  +  r  ^M^(X). 

k  X-l  3  3  Jf=l 
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Similar  considerations  to  those  above  yield 


(2.19) 


lim 
k->  <= 


Var  Y  (k). 
a  A 

Var  Y,  (k)^ 
b 


T)  >  0. 


For  the  random  car  size  model  in  one  dimension,  define 

(2.20)  Z(x)  =  total  number  of  cars  of  length  D,  where  D  has 
distribution  F,  density  f,  which  may  be  parked 
on  a  curb  of  length  x  . 


Define 

(2.21) 


M(x)  =  E  Z (x) 


1* (s )  =  f  e  SXM(x)  dx 
J0 
CO 

cp(s)  =  f  e  £Xf  (x)  dx. 

*•'  0 
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It  may  be  shown  with  the  prior  model  that  for  s  1  0,  a  v,  X 


may  be  computed  which  satisfy,  for  s  J  0, 


(2.22) 


L(s) 


-2 

vs 


d-  Xs 


Define,  for  x  >  0 

(2.23)  X(x)  =  vx  +  X 

m(x)  =  =  M(x)  -  X(x) 

M(2)(x)  =  E(Z(x)-X(x))2 

By  the  methods  of  earlier  paragraphs,  it  follows  that, 
for  a  >  0, 


(2.24) 


M(2)(a)  =  £[%(*)  dx  j"V2)(t)  dt 


+ 


ra  ra-x 

f(x)  dx  :  m(t)m(a-x-t) 

J0  J0 


If  E(D2)  <  ®, 
(2.25) 


it  follows  that 


lim  x  _1M^(x)  =  p  >  0. 

X-»  co 
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3 .  Higher  Moments. 


Let 


(3.1)  X(x)  =  total  number  of  unit  length  segments  (cars)  which  may 
be  parked  on  [0,x)  in  accord  with  a  Renyi  model,  where  each  unit  segment 
is  placed  uniformly  at  random  and  parked  only  if  it  does  not  overlap  a 
previously  parked  car.  See  (Solomon  (1966)). 

It  may  be  shown  (Solomon  (1966))  that  for  x-»co 

(3.2)  EX(x)  =  m(x)  JL  (x)  =  cx  +  c  -  1 

where  c  .748. 

Define,  for  n  >  1 


(3.3) 

M(n)(x)  =  E (X (x)  -£(*)) 

and 

M(1)(x)  s  m(x) 

It  is  known  that  (Solomon  (1966)) 

(3.4)  lim  x  Hl^\x)  =  lim  x  Krar  X  (x)  =  a  >  0. 

X  -»  eo  X  ->  ® 


Define 


(3.5)  X ' (x) ,  X" (x)  to  be  I.I.D.  as  X(x). 
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Theorem  1. 


(3.6) 


lim  x  nH^n^  (x) 


_  ,(2*1)1  a2n 

on  i 

/  n . 


and 


lim  x-(2n+1)/2M(2rH-1)(x)  =  0 


Proof.  Conditional  on  the  placement  of  the  first  unit  car  on 
(t,t+l),  x>t+l>t>0,  considerations  of  (Weiner  (1978))  yield 

that  for  asymptotic  results, 


(3.7)  E  (X  (x+1)  -  X  (x+l))n  «  E  ( [X ' (t)  -X(t)]  +  tX"(x-t)  -X(x-t)])n, 

where  t  is  chosen  uniformly  on  [0,x],  Expanding  the  right  side  by  the 
binomial  expansion  and  taking  expectations  for  n  =  3  yields 


(3.8)  M(3)(x+1)  =  6x_1  f  M(2)(t)m(x-t)  +  2  I*  H(3)  (t)  dt 

J  n  x  J  n 


and  taking  Laplace  transforms,  with 


(3.9) 


L^n^(s)  =  eS  F  e  (x)  dx, 

w0 


R 


(n) 


(s)  =  J  e  SXM^(x)  dx, 

v  r\ 


it  follows  that 


(3.10) 


(L(3)(s))'  +  L(3)(s)  =  -6R(2)(s)R(1)(s) 

s 
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with  solution 


(3.11)  l/3\s)  =  e>;p(^2J  du^)lJ  exp(-2jr  ^7"  du^jf 6R^2  ^  (v)R^ ^ (v)  dvj 


00  -u 


CD  ->U 


Since  for  0  <  x  <  1, 


(3. 12) 


CD  -n 

r  r  e  ”\ 

eXpV2J  ~7~  duy  =  _2Y  -  2  in  s  +  g(s), 
s 


where  y  =  Euler's  constant  and  g  (s)  -*  0  as  s  ■-*  0, 


(3.13) 


lim  s2R(2)(s)  =  a2 
s  1  0 


lim  s2R^^  (s)  =  0, 
s  i  0 


it  follows  that 


(3.14)  lim  s2L(3)(s) 
s  1  0 


_  co  — u 

=  J  j^exp(^-2j  du-2y)  T6r^2 ^  (v)R^  (v) Jdv  =  constant* 


For  n  =  4, 


(3.15)  M^'(x+1)  =  (J)x'1  JXH(“/(t)M^^(x-t)  dt  +  q(x) 


,(2),.w(2), 


0 


+  *  rV4)(t)  dt 

"V  1  '  ' 

x 


where 


(3.16)  q(x)  is  of  lower  order  of  x  than  x  1  F  M^2'*  (t)M^2^  (x-t)  dt 

J0 

by  the  results  for  n=l,2,3. 
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(A) 

Again  taking  Laplace  transforms,  the  resulting  equation  for  L  (s) 
has  solution 

(3.17,  l<4>(s)  -  r  (if  V  au)Tf[e*p(-2f  V  <R(2)(v)  >2+t  (v> 

s  S  V 

where  t(v)  is  the  Laplace  transform  of  q(x). 

Since 

(3.18)  lim  v2r(2)(v)  =  a2 

v  1  0 

A 

lim  v  t(v)  =  0, 
v  1  0 


if  one  denotes 

i 

(3.19)  f(s)  =  J  [exp(-2j  ^  du)j6(R(2)(v))2  +  t(v)j  dv, 

sL  V 


it  follows  by  (3.18)  and  L'Hospital's  rule  that 


(3.  20) 


lim  sf(s) 
s  1  0 


s  =  0 
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4  2v 
a  e 


Ey  (3.18),  (3.20) 


(3.21)  lim  s^L^^(s)  =  6a*. 

s  1  0 

(4) 

Since  M  (x)  is  expressible  as  a  sum  of  terms  each  increasing  in  absolute 
value,  a  Tauberian  theoren  (Kidder  (1946))  yields 

(3.22)  lim  x  2M^^(x)  =  3a*. 

x  —  a> 


dv 
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Assume  the  result  of  the  theorem  for  r=l ,  2 ,  . .  . ,  211^  moments  by 
induction.  By  the  binomial  expansion, 


(3.23)  M(2n+1)(x+l)  =  -  [*  £  (2r^1)M^')(t)M(2&fl"£) 

x  J0  1=1  £ 


(x-t)  dt 


2  T  w(2rri-l) 
x  j 


+  -  !  M'  '  (t)  dt. 


0 


(1  \ 

In  the  first  integral  on  the  right  side,  if  Mv  (t)  is  an  even  moment, 
and  the  induction  hypothesis  yields  that 


(3.24)  M(£)(t)M(2n+1“£)(x-t)  <  0(t'£/2(x-t)2n“£/2)  <0(xn). 


and  when  Laplace  transforms  are  taken  in  (3.19),  then  a  simple  computation 
by  L'Hospital's  rule,  and  an  Abelian  and  Tauberian  theorem  using  (3.20) 
yields  that 


(3.2  5) 


M(2ahl)(x)  <  o (xn ) . 


This  suffices  for  the  odd  moments. 

For  the  (2n+2)t^1  moment,  the  binomial  expansion  yields 


(3.26)  M(2n+2)(x+l)  =  -  fX  (2?t2)M(2'e)(t)M(2r>f2~2£)(x-t)  dt 

XJ0H 


2  r  (2n+2)  ,  .  _ 

+  -  ;  (t)  dt  +  q(x), 

x  „  Q 

where  q(x)  is  the  sum  of  terms  which  are  multiples  of  convolutions  of  odd  moments. 


>nd  are  of  lower  order  than  the  first  term  on  the  right  side  of  (3.22)  by  induction 
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Multiplying  by  x  and  taking  Laplace  transforms  one  obtains  an 
equation  with  solution 


(3.27)  L(2l^2)(s)  «  [exp^J  du)]  J  [exp(-2j  du)] 

S  S  V 

X  [  2  (2^2)R(2£)(v)R(2rH’2'2£)(v)  +  t(v)]  dv 


CO  -u 


where  t(v)  is  the  Laplace  transform  of  q(x). 
Let 


=>  -u 


(3.28)  f(s)  =  J  exp(-2j  pp  dupj  £  (2J^2)R(2£)  (v)R(2lH_2"2X)  (v)  +  t(v)  dv. 


-Jt=  1 


Then  by  (3. 12),  L'Hospital's  rule,  and  the  induction  hypothesis. 


(3.29)  lim  snf (x)  =  lim 

s  1  0  s  1  0  (— ) 

n 

s 


1  y  .2n+2  (2L)ia21  #  (2m-2-2£)  \a2n+2~2i  2y 

n  2SL  '  *  2iH-1-X 


Hence,  by  (3.12),  (3.29), 

(3.30)  lim  si^2l(2iH-2)(s) 
s  1  0 


By  a  Tauberian  theorem. 


1  y  2nf2  (2&)\a2£  (2n+2-2£) la2n+2~2£ 

D  A  2X  2*  *  2Itfl-X 


(3.31)  lim  x"(n+1)M(2tH'2)(x) 

X  -•  * 


1 

(n+1)  In 


£  ( 

x=i 


2n+2 
2  L  ' 


2  ( 

(2X1  lor  (2n+2-2£) 

X  _n+l-i 
2  *2 


,  2n+2-2X 

.a 


(2rrt-2)  lor 


2n+2 


2&+1(n+l)I 


and  the  theorem  is  proved. 

It  is  clear  by  the  nature  of  the  computations  that  the  moment  result 
holds  for  the  three  models  with  either  the  Renyi  or  Solomon  parking 
mechanism. 

Theorem  2.  For  any  one-dimensional  model  considered,  for  either 
the  Renyi  or  Solomon  parking  mechanisms,  where  c,  a  are  specific  to 
each  case, 

(3.22)  lim  eT — -T)  lim  l|  era^-nl— =  e"u  12  . 

L  ^  Vv-ar  *00  V  *  —  L  V  V 

2 

Proof .  The  limiting  moments  of  Theorem  1  are  those  of  a  N(0,a  ) 
random  variable,  which  uniquely  determine  the  normal  law. 


4.  Two  Dimensions. 

Let 

(4.1)  R(x,y)  =  total  number  of  unit  squares  which  may  be  parked  on  an 
xx  y  rectangle  in  the  Renyi  model  in  two  dimensions. 

Lemma.  For  x,y-»co,  all  n  >  1 


(4.2)  (i)  E(X(x)X(y)  -X(x)£(y))2n  ~  f  (ca)2n(xy)n(x+  y)n 


nl  2 


2n+1  n+2 

(4.2)  (ii)  E(X(x)Y(y)  -i(x)X(y))  ~  o((xy)  (x+  y))  ' 


with  Y(y)  I.I.D.  as  X(y). 
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Proof. 


For  x,y  large,  all  n  , 


(*•3)  E[Y(y)(X(x)  -X(x))  +  X(x)(Y(y)  -X(y))]n. 

The  independence  of  X(x),Y(y),  the  asymptotic  expressions 

(4.4)  (i)  Si  (x)  ~  cx 

(4.4)  (ii)  E(Y(y))n  =  E[(Y(y)  -X(y))+X(y)]n  ~  Ce(y))n~cnyn. 

a  binomial  expansion  of  (4.3),  Theorem  1,  an  explicit  asymptotic 

computation  of  (4.2) (i)  for  n=l,2,  (4.2)(ii)  for  n=l,  and  an  induction 
establish  the  lemma. 


(4.5) 


lim  E 
x,y-»oo 


{ 


exp 


ur 


R(x,y)  -  a  (x)i (y) 


'E(R(x,y)  -£(xU(y))i 


*  llm  E{exp[~u( 7^=  — )]}  ' 


-u2/2 


x,y-*» 


2  2 

c  a  xy(x+y) 


Proof.  For  x,y-4»,  n  >  1, 

(4.6)  E(R(x,y)  -X(x)X(y))n  =  E[(Y(y)(X(x)  -X(x))  +  (X  (x)  (Y (y)  -X(y)))]Q. 
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The  lemma  immediately  yields  that  for  x,y-»a>,  and  (4.6),  that 


(4.7) (i) 


E (R(x ,y)  -  l (x)j (y ) ) 2n  ^  (2n)I 
(ca)2n[xy(x  + y)]n  2nn! 


(4.7)(ii) 


E(R(x,y)  -Si  (x)j  (y)) 


2n+  1 


[xy(x+  y)] 


n  + 


0, 


all  n  >  0. 


The  right  sides  of  (4.7)  (i),  (4.7)  (ii)  are  the  2nt^1  and  (2n+ l)St 
moments  of  a  N(0,1)  random  variable,  which  determine  the  distribution, 
proving  the  theorem. 

Similar  results  hold  for  the  Solomon  model  and  other  models  con¬ 
sidered  earlier.  A  bivariate  central  limit  theorem  holds  for  alternating 
squares  or  rectangles  models. 
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